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In this paper we derive the quantum statistical and dynamical properties of non- 
linear optical couplers composed of two nonlinear waveguides operating by the sec- 

Qh' ond subharmonic generation, which are coupled linearly through evanescent waves 

-t— > ■ 

and nonlinearly through nondegenerate optical parametric interaction. Main at- 
d | 

tention is paid to generation and transmission of nonclassical light, based on a 
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discussion of squeezing phenomenon, normalized second-order correlation function, 



£> ' and quasiprobability distribution functions. Initially coherent, number and thermal 

(N 

states of optical beams are considered. In particular, results are discussed in depon- 
es 

dence on the strength of the nonlinear coupling relatively to the linear coupling. We 

show that if the Fock state |1) enters the first waveguide and the vacuum state |0) 



enters the second waveguide, the coupler can serve as a generator of squeezed vacuum 
state gevorned by the coupler parameters. Further, if thermal fields enter initially 



the waveguides the coupler plays similar role as a microwave Josephson-j unction 
parametric amplifier to generate squeezed thermal light. 
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I. INTRODUCTION 

In quantum optics many simple quantum systems have been examined from the point 
of view of completely quantum statistical description including not only amplitude and 
intensity (energy) development of such systems, but also higher-order moments and complete 
statistical behaviour. Such results have fundamental physical meaning for interpretation of 
quantum theory [1] and they are useful for applications in optoelectronics and photonics as 
well. These results can be successfully transferred to more complicated and more practical 
systems, such as optical couplers composed of two or more waveguides connected linearly 
by means of evanescent waves. The waveguides used can be linear or nonlinear employing 
various nonlinear optical processes, such as optical parametric processes, Kerr effect, Raman 
or Brillouin scattering, etc. Such devices play important role in optics, optoelectronics and 
photonics as switching and memory elements for all-optical devices (optical processors and 
computers). When one linear and the other nonlinear waveguides are employed, we have a 
nonlinear optical coupler producing nonclassical light in the nonlinear waveguide which can 
be controlled from the linear waveguide, i.e. one can control light by light. The generation 
and transmission of nonclassical light exhibiting squeezed vacuum fluctuations and/or sub- 
Poissonian photon statistics in nonlinear optical couplers can further be supported when 
all the waveguides are nonlinear. The possibility to generate and to transmit effectively 
nonclassical light in this way is interesting especially in optical communication and high- 
precision measurements where the reduction of quantum noise increases the precision. In 
the present paper nonlinear couplers have been examined composed of linear and nonlinear 
waveguides [2] (and references therein) with a particular attention to quantum statistical 
properties of such devices |3ru| related to quantum noise properties. These devices are 
useful for generation and transmission of nonclassical light and new interesting effects can 



performed 
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4lQ|. Also Schrodinger-cat states can be 



be obtained if phase mismatches are involved 

transmitted through nonlinear couplers [U| and stability analysis of such devices can be 



Nonlinear codirectional and contradirectional couplers composed of two nonlinear waveg- 
uides operating by second harmonic generation or bynondegenerate optical parametric pro- 
cesses can exhibit interesting switching properties [l3|, [14| . Quantum-consistent description 
of contrapropagating beams can be developed, which permits to formulate the problem in 



the Hamilton formalism Il5|. Phase mismatches inside the nonlinear waveguides and be- 
tween them can be taken into account |10(]. Interesting results can be obtained for the 
quantum statistical properties of nonlinear optical couplers operating by means of Raman 
and Brillouin scattering |3|. 

In this paper we continue in investigation of quantum statistical properties of nonlin- 
ear couplers composed of two waveguides operating by the second subharmonic generation 
assuming strong coherent pumping and linear exchange of energy between waveguides by 
means of evanescent waves, however we additionally take into account the influence of non- 
linear coupling of the parametric type of both the waveguides. In section 2 we describe 
dynamics of the system under discussion together with the solution of the equations of mo- 
tion, in section 3 we derive squeezing characteristics of generated light, section 4 is devoted 
to a discussion of sub-Poissonian statistics, section 5 includes results for quasidistribution 
functions and finally we summarize main conclusions in section 6. 

II. MODEL DESCRIPTION AND EXACT SOLUTION 
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Let us consider a system described by the Hamiltonian H such that, 

2 
\_] \ ^jO-j^j + -\? 0*1 BXp(i/jLjt) + h.C. > 

3=1 

+A 3 < hia\ exp[i(f>i(t)] + h.c. > + A 4 {aid 2 exp[— ifaif)] + h.c.} , (1) 

where a\ (oj), &2 (^2) are annihilation (creation) operators of the fundamental modes in 
the first and second waveguides having frequency U\ and 0J2, respectively, (j,j are related 
with the frequency of the second-harmonic modes described classically as strong coherent 
fields, 4>j(t), j = 1,2, are related to the difference- and sum-frequencies of modes 1 and 2, 
respectively, Ai and A2 are nonlinear coupling constants for the second subharmonic gener- 
ation in the first and second waveguides, respectively, A3 is the coupling constant for linear 
exchange between waveguides through evanescent waves, A4 is the coupling constant for the 
nonlinear exchange through simultaneous annihilation or creation of a photon in both the 
subharmonic modes on expense of pumping and h.c. means the Hermitian conjugate terms 
(for further details concerning the optical parametric processes, see [l6( (Chap. 10)). When 
fij = and only the degenerate term is considered, we have the well-known Hamiltonian, 
in the interaction picture, for squeezed light generation [17| . where Ai (or A 2 ) represents 



the coupling constant proportional to the quadratic susceptibility, of the second-order non- 
linear process (degenerate parametric down-conversion with classical coherent pumping), or 
the coupling constant proportional to the cubic susceptibility, of the third-order nonlinear 
process (degenerate four- wave mixing with classical coherent pumping) [18]. If additionally 
(j>i(t) = 4>2{t) = 0, the Hamiltonian (1) represents a mixture of second subharmonic genera- 
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tion, frequency conversion and parametric amplification in the interaction picture 

It is important to mention that we treat the problem of propagation in the Hamiltonian 
formalism neglecting dispersion. Thus if case all waves are propagating with the same 
velocity, time t and space z relate by the velocity of propagation v, z = vt. Schematically, 
this Hamiltonian is represented in Fig. 1. 



a 2 (o) 



ai(0) 





z = vt 



Fig.l Scheme of quantum nonlinear coupler with linear and nonlinear coupling formed 
from two nonlinear waveguides described by the quadratic susceptibility x^ ■ The beams 
are described by the photon annihilation operators as indecated; z = vt is the interaction 
length. Both the waveguides are pumped by strong classical coherent waves. Outgoing fields 
are examined as single or compound modes by means of homodyne detection to observe 
squeezing of vacuum fluctuations, or by means a set of photodetectors to measure photon 
antibunchibng and sub-Poissonian photon statistics in the standard ways. 

In fact the Hamiltonian (1) can be regarded as a generalization of the models given in refs. 
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24j. For example, if we take both Ai and A2 to be zeros, then we shall be left with the 



Hamiltonian which describes the back-action evading amplifiers, where the Hamiltonian in 
this case can be constructed by combining parametric amplifiers and parametric frequency 
converters with two different coupling parameters. On the other hand, if we take fij = 
and drop the time dependent phases, then the Hamiltonian (1) will be consistent with the 



Hamiltonian given in ref. 25[, where the wave functions for both the number state and 
coherent state and the Green's function have been obtained. It is also interesting to point 
out that the Hamiltonian (1) contains ten generators based on the group sp(4,R), which 
represents the most general type of the two-mode quadratic Hamiltonian [26|. This will 
enable us to reconsider the problem from Lie algebra point of view, where the most general 
solution for the wave functions may be obtained. For more details, see for example refs. 
27H29|. where the wave function for some special cases of the above Hamiltonian has been 
obtained using the Lie algebra technique. 

Annihilation and creation operators satisfy the boson commutation relations 



[ai,a/] 



-»!tfj 



where 0^ is the Kronecker delta. 



The equations of the motion in the Heisenberg picture for the Hamiltonian (1) are 
— — = —iuiai — 2i\ia\exp(itfii) — 1X30,2 exp[—icj)i(t)] — iA4a 2 exp[z0 2 (£)], 

— — = — io> 2 a 2 — 2iX2al exp (it /j, 2 ) — iX^aiexp[i(j)i(t)] — zA 4 a{ exp[i(/> 2 (£)]. 



(2) 

(3a) 
(36) 



Substituting a\ = Aexp(~/ii) and a 2 = -Bexp(^/i 2 ), slowly varying forms of the operators, 

having the operators Oj as well as A and B time dependent, equations (3) take the form 

dA ., /ii , 



dt 



-i(u)i + —)A - 2iX 1 A ] - iX 3 Bexp 



-iX A & 



cxp 



.{Hi + fl>2)t 



'.(fi2-til)t . , , ' 

% l(j>l{t) 



i<fo(t) 



dB /i 2 ,- 



iX^As exp 



2iX 2 B ] -iX 3 Aexp 

-l h l<j>2\t) 



i ^-^ +iMt) 



(4a) 



(4b) 



Equations (4) with their Hermitian conjugates give a close system of four differential equa- 
tions with time-dependent coefficients which cannot be solved directly and hence some 
restrictions should be considered, so that we shall consider 4>i(t) = |(/i 2 — fj,i)t and 
4>2(t) = |(/^2 + Hi)~t- Then the solutions of the system (4), which yield the relations be- 
tween input and output modes, can be obtained, after some tedious calculations, as 

Oi(t) exp(-ztf ) = ai(0)ifi(£) + oI(0)Li(t) + a 2 (0)M 1 (t) + c&(0)Ni(t), (5a) 

a 2 (t) exp(-itf) = a 2 (0)K 2 (t) + a 2 (0)L 2 (t) + oi(0)M 2 (t) + al(0)A^ 2 (t), (56) 

where the time dependent coefficients, which contain all the features of the structure, are 
given by 



K 1 (t)=F 1 (t)- 



[A; + + A;_]G' 1 (t) + [A + ^ + A_] 1 S(i) 
°i 



(6a) 



Li(t) = -- 



Mi(t) = - 
Ni(t)^- 



[k+ - fc_]Gi(t) + [\ + - - \.}S(t) 
9i . 
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whereas 



- ^ ) C(t) 
9i J 

g A en) 

9i J 
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— i 



92 



K 2 (t) = F 2 (t) - - 



[A + + A_]G 1 (£) + [J+^ + J_]S(t) 

[A+ - \-\Gi{t) + [J+- - J-]S(t) 
9i . 



92 



1 



[J+ + J-]G 2 (t) + [A+^ + \.}S(t) 
9i 

92, 



[J + -J.]G 2 (t) + [\ + -\-^]S(t) 

9i 



M 2 (t) = - 



N 2 (t) 
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l+>j±)C(t)-i 
In the above equations we have defined 



[\+ + \-}G 2 (t) + [k + + k-^]S(t) 

9i 



1 C(t) - % 



[A+ - \-}G 2 (t) + [k, 



92- 

9i 
92-, 



k_^]S(t) 
9i 



A± — A3 ± A4, 

k± = u)! + -fii ±2Ai, 

J± = u 2 + -/i 2 ± 2A 2 , 

9l = k_\ + + X-J+, 

g 2 = k + \- + X+J-, 

and 

F t (t) = cos(ifii) cos 2 9 + cos(t0 2 ) sin 2 9, 

F 2 (t) = cos(tQ 2 ) cos 2 9 + cos(tf2i) sin 2 9, 
sin(tfii) 2 sin(tfi 2 ) . 9 „ 

Gi(t) = =; COS 6M 

\l\_ J4 2 _ 

""'"' sin 2 0, 



sin(tfl 2 ) 2 
G2(t) = — = cos 9 



Qo 



n 2 

sin(tf2i" 



C(t) = -.r^\cos(tn 2 ) - cos(iQi)l sin(20) 
2 V #2 

sin(£f2 2 ) sin(tfii 



s ^m 



sin(26>), 



fi 2 fii 

where we have introduced the abbreviations 
-1 / ^V9i92 



tan 



J_J, 



fv /Cj 



n 2 



[fi 2 cos 2 + fi 2 sin 2 - v /#7J 2 ~sm(20)] ' 
[fi| cos 2 + fi 2 sin 2 6> + ^/gTg^ sin(29)] ' 



with Qf = A_A + + k_k + and fi 2 . = A_A + + J_J^ 



(6b) 
(6c) 

(6d) 

(7a) 
(7b) 
(7c) 
(Id) 



(8) 

(9a) 
(9o) 
(9c) 

(9d) 

(9e) 

(9/) 

(10a) 

(10b) 
(10c) 



One can see from this solution that when fii and fi 2 are real, the coupler switches the en- 
ergy between the modes which propagate inside since the solution will include trigonometric 



functions [19|. Nevertheless, if f^ and f2 2 are pure imaginary, the Heisenberg solutions at- 
tribute hyperbolic functions, which are growing rapidly with time, and the coupler operates 
as amplifier for the input modes 30[. So that the behaviour of the coupler will be indicated 
essentially by the relation between coupling constants. 

For the time dependent coefficients, we can easily obtain the following relations 



|^.(t)| 2 + |M,(t)| 2 = l + |L,(t)| 2 + |ty(t)| 



2. 



3 



1,2, 



(11a) 
(116) 

(lie) 



K x {t)N 2 {t) + M x (t)L 2 (t) = N^K^t) + L x (t)M 2 (t), 

Ki(t)M 2 *(t) + M 1 (t)K;(t) = Li(t)iVJ(t) + N^L^t), 
in correspondence to boson commutation rules (2). 

In what follows, we shall employ the results obtained in the present section to treat the 
squeezing phenomena, normalized second-order correlation function, as well as quasiproba- 
bility distribution functions for the model under consideration. 



III. SQUEEZING PHENOMENON 



Squeezing is a pure nonclassical phenomenon and squeezed states have less noise in one 
field quadrature than a coherent state. On the other hand, this means that there is an excess 
of noise in the conjugate quadrature, since the product of canonically conjugate variances 
must satisfy the uncertainty relation. This light has a lot of applications, e.g. in optical 
communication networks [31[, in interferometric techniques [321 ] . and in optical waveguide 
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33j. Generation of squeezed light has been observed in many optical processes [34|, 135 ]. 



Investigation of the squeezing properties of the radiation field is a central topic in quantum 

optics which can be measured by homodyne detection where the signal is superimposed on 

a strong coherent beam of the local oscillator. 

For this purpose we define the position and momentum operators, which are related to 

the conjugate electric and magnetic field operators E and H of electromagnetic field, for 

each mode in terms of a,j(t) and a At) as 

, itfXj 
Y 



Mt) = o 



YAt) = — 



«i(*) ex P(^) + a](t) exp( 



—itpij 



a,j(t) exp( — -) — dj\t) exp 



2 ' 
—itfij 



(12) 



(13) 



2 ' 3 w ~"" v 2 
where we have considered Hj(t) to be the phase of the local oscillator, without loss of 

generality, to cancel the high frequency terms, and j = 1,2 stands for mode 1 and mode 2, 



with AXj(t) 
One of the fo 
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respectively. These operators satisfy the commutation relations 

'x j (t),Y j (tj\=±, (14) 

so that the uncertainty relations are 

AXj(t)AYj(t) > \, (15) 

i 

((A^(t)) 2 )] 2 = [(X|(t))-(l,(t)) 2 ]i 
lowing squeezing conditions for each mode can occur, 
S 3 (t)=A((AX 3 (t)f)-l<Q, 

Q,(t)=4((Ay,(t)) 2 )-i<0, (16) 

i.e. negative values of these quantities express squeezing of vacuum fluctuations. Here we 
study squeezing phenomenon when the modes are initially prepared in thermal-states (or in 
number states since both of these two cases, number states and thermal-states, have identical 
quadrature variances) with the average thermal photon numbers fij, j = 1,2 as well as 
in the coherent states. More details on the evolution of thermal light in the model under 
discussion will be adopted in section 5. Now for the quantities Sj(t) and Qj(t), provided that 
both the modes are initially in the thermal states, we have for the first mode the following 
expressions 

S^t) = 2n 1 [|L 1 (t)| 2 + l^l 2 ] + 2n 2 [\N 1 (t)\ 2 + |Mx(t)| 2 ] + 2|L 1 (£)| 2 + 2|iV 1 (t)| 2 
+2(2ni + l)[Li(£)^i(£) + c.c] + 2(2n 2 + l)[Mi(t)iVi(t) + c.c], (17) 

Q 1 (t) = 2n 1 [\L 1 (t)\ 2 + |^(£)| 2 ] + 2n 2 [|iV 1 (t)| 2 + |Mx(t)| 2 ] + 2|L x (t)| 2 + 2|iV 1 (t)| 2 

-2(2^ + l)[Li(t)Jri(t) + c.c] - 2(2n 2 + l)[Afi(i)Wi(t) + c.c], (18) 

where c.c. means the complex conjugate terms. The corresponding expressions for the 
second mode can be obtained from (17) and (18) by using the interchange 1 -H- 2. However, 
the other expressions related to the injected coherent light initially in the coupler are the 
same (17) and (18) but just put fij = 0. 

It is known that the nonlinear coupler is a source of optical fields, the statistical prop- 
erties of which are changed as a result of the linear and nonlinear interaction inside and 
between waveguides. Consequently, one can generate nonclassical light from one input and, 
in addation, it can be switched. 

We have plotted Si(t), Q\{t) in Figs. 2a,b and S 2 (t), <5 2 (^) i n Figs. 3a,b, when the initial 
light is coherent, for different values of A&. Further we have chosen A3 = 1 for all curves 
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FIG. 2: Squeezing phenomenon for mode 1 when the modes are initially in coherent light and 
in thermal light. For initial coherent light : a) for the first component S\(t); b) for the second 
component Qi(t); A3 = 1 for all curves; curves A, B and C are corresponding to Ai = A2 = A4 = 
0.25, Ai = A2 = A4 = 0.2 and Ai = 0.17, A2 = A4 = 0.2, respectively. For initial thermal light: 
c) the first component S\(t) with n\ = 0.5, n,2 = 0.5 (solid curve), 1.5 (dashed curve) and the 
coupling constants Xj are the same as those for the curve C when the light is initially coherent; 
straight line has been put to show the bound of squeezing. 
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and for the curve A : Ai = A 2 = A 4 = 0.25; for the curve B : Ai = A 2 = A 4 = 0.20, 
and for the curve C : X± — 0.17, A 2 = A 4 = 0.2. On the other hand, Fig. 2c gives Si(t) 
(first mode) when the initial light is thermal light with coupling constants as those for the 
curve C, where fi\ = 0.5 and n 2 = 0.5 (solid curve), 1.5 (dashed curve); and straight line 
shows the bound of squeezing of the curves. Firstly, we start our discussion by studying the 
case of input coherent light. From these figures we can see how the coherent states, which 
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FIG. 3: Squeezing phenomenon for mode 2: a) for the first component ^(i); b) for the second 
component Q 2 (t); the values of the parameters A& are as in Fig. 2. 

are minimum- uncertainty states, evolve in the coupler to produce squeezed light. We can 
observe the oscillatory behaviour in these curves, showing that squeezing can be switched 
from one waveguide to the other in the course of time during power transfer. Moreover, 
squeezing can be interchanged between the two quadratures of the same waveguide. More 
precisely, for mode 1, squeezing can occur for all selected values of A^ in Si(t), but in 
Qi(t) only curves A, B can exhibit squeezing, as shown in Figs. 2a,b, which reflects the 
dependence of nonclassical behaviour on the strength of subharmonic generation. For mode 
2 we can see squeezing in all curves in both the quadratures, as shown in Figs. 3a,b. It 
can be easily seen that the amount of squeezing is sensitive to the strength of coupling A& 
and that in general its values in the second component are more pronounced than those in 
the first one. Now if we turn our attention to the case of injected thermal light, i.e. Fig. 
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2c, we can observe that squeezing is available in the large interaction time. Further, Si(t) 
exhibits oscillatory behaviour and it evolved from unsqueezed values in the short range 
of interaction time, owing to the fact that thermal-states are not minimum-uncertainty 
states, into squeezed values and eventually unsqueezed values can be recovered. Indeed, we 
noted numerically that this behaviour is periodically recovered with the time. Moreover, by 
comparing the dashed curve with the solid one, we can see that increasing of the photon 
number in the second waveguide causes decreasing of the amount of squeezing in the first 
one. This is related with the effect of evanescent waves bewteen waveguides and shows how 
one can control light by light in the coupler. Finally, we can conclude that by controlling 
the input average thermal photon number and the interaction time (or on the length of 
the coupler), the interaction under consideration can generate squeezed thermal light. It 



is worthwhile to refere to 



where more discussions related to squeezed thermal states 



are given. Furthermore, squeezing of thermal radiation field has been already produced 



m 



where a thermal input field 



in a microwave Josephson-j unction parametric amplifier 

has been introduced to the squeezing device and the generated field has exhibited noise 

reduction. 

IV. SECOND-ORDER CORRELATION FUNCTION 

Starting with the experiment of Hanbury Brown and Twiss, strong interest in the photon- 
counting statistics of optical fields began. Traditional diffraction and interference experi- 
ments and spectral measurements may be considered as being performed in the domain of 
one photon or linear optics. The theory of higher-order optical phenomena, describe d by 
higher-order correlation functions of the electromagnetic field, was founded by Glauber 38] , 
who introduced the measure of super-Poissonian statistics (classical phenomenon) and sub- 
Poissonian statistics (nonclassical phenomenon) of photons in any state, which is given by 
the normalized normal second-order correlation function defined as 

(2 )(t) = (*?(t)tij(t)) 

i ((An J -(t)) 2 )-(qt(t)a J -(t)) 

(a]m(t)) 2 ' ( , 

where the subscript j relates to the jth mode and ((Afij(t)) 2 ) are the photon number 
variances, which can be obtained from the relation 
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<(An,(t)) 2 > = <(at(t)a,(t)) 2 > - (a](t)%(t)} 2 . (20) 

Then it holds that g, (t) < 1 for sub-Poissonian distribution of photons, g): (t) > 1 for 
super-Poissonian distribution of photons and when g- (t) = 1 Poissonian distribution occurs. 

(2) 

The degree of coherence g, (£) can be measured by a set of two detectors. An application 
of radiation exhibiting the sub-Poissonian statistics to optical communications has been 



considered in [39| 



The most familiar quantum states from the earlier days of quantum mechanics are co- 
herent and number states. Following the development of the quantum theory of radiation 
and with the advent of the laser, the coherent states of the field, that mostly describe a 
classical electromagnetic field, were widely studied. These states are minimum-uncertainty 
states and have Poissonian distribution of photons and they may be evolved in the non- 
linear optical coupler to generate nonclassical light. On the other hand, number states are 
purely nonclassical states (they always exhibit sub-Poissonian statistics) and there is great 



interest for their preparation and quantum non- demolition detection 40h42| . because they 
exhibit the maximum channel capacity, i.e. they provide the maximum of information that 
can be transmitted by a single photon, and the minimum time-energy product in optical 
communications 43]. 



Here we shall study the intensities of the fields as well as the normalized normal second- 
order correlation function for mode 1 when both the modes are initially in the coherent 
states |a)i, \/3)2 or in the number states |n)i, 1 771)2- Then the photon number variance in the 
coherent state is given by 

((An,(t)) 2 ) coh = [V^ )2 {t) +4|V^(t)| 2 ]M 2 + [Kp(t) +4|V^(t)| 2 p| 2 

+[\V 7 U \t)\ 2 + \V®{t)\*\(\°\ 2 + W\ 2 ) + [\V 6 (J \t)\ 2 + 2|^ ) (i)| 2 + 2|K^(t)| 2 ] 
+ {a 2 [2V 1 u \t)V; U \t) + V; U \t)V^\t)} + f3 2 [2V 2 U \t)V* iJ \t) + V 7 u \t)V 6 U \t)] 
+ap[vi J) (t)V^ j) (t) + 2V; U \t)V 7 U \t) + 2V 5 * U \t)V 7 * U \t) + V 2 U) (t)V 6 U) (t)] 
+a*P[vi J \t)vi' J \t) + 2V} j \t)V 6 (j \t) + 2V* (j \t)V£ {j \t) + V 2 U) (t)V 7 iJ \t)} + c.c.},(21) 
while the expectation value of the photon number is 
(a](t)a,(t)) coh = \a\ 2 V{ j \t) + |/3|V 2 °' } (t) + V^\t) 

+ \a 2 V; U \t) + f3 2 V* U \t) + a*(3V 7 {j \t) + af3V 6 U \t) + c.c.j . (22) 

For initial number state we find the photon number variance in the form 
((An,(t)) 2 ) n = 2|I/ 4 (J) (£)|V + n + 1) + 2\V 5 U \t)\ 2 (m 2 + m + 1) 

+(\V^(t)\ 2 + \V 7 U \t)\ 2 )(n + m + 2mn), (23) 
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while the expectation value of the photon number is 



it/ 



-U) 



(a (*)&,-(*)>„ = n^ w (t) + mV^{t) + V^>(t) 



U), 



rU), 



where 



v} 

v! 

v} 



(t) 
(t) 
(t) 
(t) 



\K 3 {t)\ 2 + \L 3 (t)\\ 

\M 3 {t)\ 2 + \NM\ 
|iV,(t)| 2 + |L,(i)| 2 , 



(t) = M;(t)N,(t), 

{t) = K 3 {t)N*{t) + L*(t)M J (t) 1 
(t) = M 3 {t)K*{t) + N*{t)L 3 {t), 



(24) 

(25a) 
(256) 
(25c) 
(25d) 
(25e) 
(25/) 
(25s) 



and j = 1,2 corresponding to first and second mode, respectively. It is important to study 
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FIG. 4: Mean photon number against time t for mode 1 (solid curve) and mode 2 (dashed curve) 
when both the modes are initially in the coherent states with a = 20,/? = 5, Ai = 0.17, A2 = 
0.2, A 3 = 1: a) A 4 = 0.2; b) A 4 = 2. 

the evolution of the mean photon numbers (intensities) (ajt)dj(t)) inside the waveguides 
of the coupler to visualize how the energy is exchanged between the waveguides. For this 
purpose we show Fig. 4 in which the mean photon number (22) of the beams is plotted 
against the time t for shown values of the parameters. The solid and dashed curves are 
related to the first and second beams, respectively. We note that the essential for the 
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behaviour of the coupler under consideration is relation of powers of the linear (A3) and 
nonlinear (A 4 ) coupling constants. To be more specific, for A 3 > A 4 (Fig. 4a), the intensities 
evelve oscillatory with time t, which means that the periodic power transfer occurs between 
waveguides and the coupler operates as an optical switcher. Further, at certain values of 
time, corresponding to intersections of the two curves, all energy in the coupler becomes 
equally shared between the propagating modes. However, for A3 < A4 (Fig. 4b), the initial 
intensities are amplified in the course of time and the coupler operates as an amplifier for 
input modes. 




(a) (b) 

(2) 

FIG. 5: Normalized normal second-order correlation function g\ (t) for mode 1 when both the 
modes are initially in the coherent states with a = 5, f3 = 20 (solid curve) and a = 20, (3 = 5 
(dashed curve): a) for both curves Ai = 0.17, A2 = 0.2, A3 = 1 and A4 = 0.2; b) for both curves 
Ai = 0.17, A = 0.2, A 3 = 1 and A 4 = 2. 



A similar behaviour is expectable for the normalized normal second-order correlation 

function for mode 1 if initially both the modes are in coherent states (Figs. 5 for shown 

(2) 
values of the parameters). In other words, for A3 > A4, we observe that g\ (t) has oscillatory 

behaviour between Poissonian and super-Poissonian statistics, i.e. coherent light can be 

approximately recovered at certain values of time. This behaviour is independent of the 

initial amplitudes of the input light (compare solid and dashed curves). On the other hand, 
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for A3 < A4, the oscillatory behaviour disappears and the fields begin to be localized in the 
waveguides into which they were initially launched. The interesting point, which could be 
realized here, is that there is a possibility to generate sub-Poissonian light from the initial 
Poissonian light input into the coupler provided that a > /3 (Fig. 5b). 
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FIG. 6: Normalized normal second-order correlation function g\ (t) for mode 1 when both the 
modes are initially in the number states with n = 5, m = 50 (solid curve) and n = 50, m = 5 
(dashed curve): a) A& have the same values as in Fig. 5a; b) A/% have the same values as in Fig. 5b. 



The situation will be quite different if we inject initially number states in the coupler, 
as is illustrated in Figs. 6, where we see that the initial sub-Poissonian statistics for Fock 
state are not recovered in the progress of time t and super- Poissonian statistics dominate. 
However, g\ (t) exhibits oscillatory behaviour under the condition provided that the linear 
coupling is stronger than the nonlinear coupling (Fig. 6a). 

We can conclude that this structure can be used to generate nonclassical light from 
classical light, e.g. coherent light, by controlling the device design and the initial input 
field. Of course, this is based on the fact that when electromagnetic fields are guided inside 
the structure, exchange of energy between the two waveguides is possible because of the 
evanescent field between the waveguides [441 ]. 
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V. QUASIPROBABILITY FUNCTIONS 



Here we shall continue in our investigation for the statistical properties of the system 
under discussion in the basis of quasiprobability distribution functions for compound modes 
when both the modes are initially in number, coherent and thermal states. 

There are three types of these functions: Wigner W-, Husimi Q- and Glauber P-functions. 
These functions give a complete description for the statistical properties of a microscopic 
system and provide insight into the nonclassical features of the radiation fields. For example, 
the density operator for the quantum mechanical system can be expressed in terms of them 
and the various moments of the system operators may be obtained by appropriate integration 



in phase space using these functions 



determined in homodyne tomography 45] 



6J. Furthermore, these quasidistributions can been 



On the other hand, as we have mentioned before, propagation of waves inside the nonlin- 
ear directional coupler causes energy exchange between the waveguides owing to the evanes- 
cent waves and hence if the measurement of an observable in the first waveguide is performed, 
this projects the state of the other waveguide into a new state; so it would be convenient 
to consider in our investigation not only the joint quasiprobability functions but also these 
functions for single modes. 

The starting point for our analysis is the s-parametrized characteristic function which 
is complex in its nature and may be used also to generate the different moments of the 
quantum system by means of differentiation. The two-mode s-parametrized characteristic 
function is given by 



C (2) (G,C 2 ,M) = Tr<Uo)exp 



(26) 



L«=i 

where s takes on values 1,0 and —1 corresponding to normally, symmetrically and antinor- 
mally ordered characteristic functions, respectively, p(0) is the initial density operator for 
the model and Tr denotes trace of the operator. 

The s-parametrized quasiprobability distribution functions are defined as the Fourier 
transform of the s-parametrized characteristic function by 

W^(a u a 2 ,s,t) = -L J J d 2 C 1 d 2 C 2 C^(CuC2,s,t)exp 

where C^ 2 \Cij C2, s, t) is given by (26). When s = 1, 0, — 1, equation (27) gives formally P- , 
W- and Q-functions, respectively. 



5^(«iC - a Ki 



1=1 



(27) 
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The corresponding single-mode s-parametrized characteristic and quasiprobability func- 
tions are 

CW(0, s,t) = Tr |p(0) exp [||C/ + &at(t) - C*%(£)] } , (28) 

W^{ aj ,a,t) = ^ J d 2 QC (1 \Q,s,t)exp( aj q -Qa*), j = 1,2. (29) 

The superscripts (1) and (2) in the above equations stand for single-mode case and two- mode 
case, respectively. 

The various moments of the bosonic operators for the system, using the characteristic 
functions and quasiprobability functions, in the normal form (N), antinormal form (A) and 
symmetrical form (S), corresponding to s = 1, —1,0, respectively, can be obtained by 

<n«rww>iw =n g cr a J ( _ g) n J g(2) (c^t)^ 1 ,_ 1 , i^^ 

r 2 

= / ^ 2 )(a, s ,t) (s=li _ li0 )n«r j «? d2 ^' ( 3 °) 

J i=i 

where n,j,mj are positive integers, C — (C11C2), « = (cti,&2), and the integral is taken 

over a.\,a2 in phase space. For example, when ni = mi = 1 and 712 = m<i = 0, then 

(a\(t)ai(t)) N = (a\(t)ai(t)), (al(t)ai(t)} A = (ai(t)a[{t)) , and (a\(t)ai(t))s 

a[(t)ai(t) + ai(t)a\(t)}. The formula (30) is valid for the single and compound modes 



l/At 



owing to the normalization of quasiprobability functions and taking into account that the 
single mode characteristic function can be obtain from that for two modes by simply setting 
one of the parameters (Ci or (2) equals zero, 
(i) Input Fock states 

It is known that the nonlinear directional coupler is an important optical device to generate 
nonclassical light in the context of control of light in the nonlinear medium. So the initial 
input light has a direct relation to the output light. In fact, investigation of output light from 
the coupler, when the number states are initially injected 43|, |46M48| . took little attention 
compared with the injected coherent states. This seems to be related to the complexity 
of calculations. However, some interesting results have been extracted by considering such 
situation 46|, |48) . For example, we can mention, in the linear directional coupler, displaced 



number states can be generated if a number state enters waveguide 1 and a strong coherent 
field enters waveguide 2 46j ; also a coherent state has been obtained in the nondegenerate 
optical parametric symmetric coupler when one of the modes enters the coupler in the Fock 
state |1) and the other modes are in vacuum states 48j. Here we shall turn our attention 
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to deduce the quasiprobability functions for the Hamiltonian (1) when the two modes are 
initially uncorrelated and enter the coupler in number states. Of course, this will give general 

formulas having wide applicability for special cases [6|, Il9l-l2l| by appropriate choice of the 

parameters. It is important to mention that some of these special cases have not been 
considered before [6|, Il9l-l21|. 

The density operator for two-mode number states is 

p n (0) = |n) 1 |m) 22 (m| 1 (n|. (31) 



Inserting (31) into (26), the two-mode s-parametrized characteristic function takes the form 

^(Kil 2 + IC 2 | 2 )-i(b(t)| 2 + Mt)| 2 ) 



C«(Ci,G«,M) 



exp 



(32) 



xL n (\ m (t)\ 2 )L m (Mt)n 
where 

77! (t) = CiKl(t) - CiLi(t) + C 2 M*(t) - <£N 2 {t), 

m (t) = (iM*(t) - QNxit) + ( 2 K*(t) - QL 2 (t), 

and L n represents the Laguerre polynomial. 

Equations (32) and (27) yield Wigner function for two-mode number states; after some 



(33a) 
(336) 



manipulations, we have the following expression 

^S(«i,« 2 ,, = 0,t) = A(_i ) (^) Ln(4 |A 1 (t)| 2 )L m (4|A 2 (t)|^ 



71^ 



,a 2l s 
x exp [-2 (lA.it)] 2 + \A 2 (t)\ 2 )] , 



(34) 



where 



Ai(t) = a t K*{t) - alL^t) + a 2 M;{t) - a*N 2 {t), (35a) 

A 2 (£) = «iM*(t) - a*N t (t) + a 2 K* 2 (t) - a* 2 L 2 (t). (356) 

Equation (34) cannot be factorized owing to the intermodal correlation between the propa- 
gating modes inside the coupler and this is clear since (34) includes terms like a\a 2 ,a\a 2) 
etc. 

The single-mode s-parametrized characteristic function for the first mode can be obtained, 
by similar way as for the two-mode case, from (28) as 



exp 



;Ki 



(h(t)| 2 + h(t)| 2 ) 



L n (h(*)| 2 )L m (h(t)| 2 ), 



(36) 



C$L(Ci,M) 

where 

v 1 (t) = ( 1 K* 1 (t)-QLi(t), (37a) 

V2 (t) = CxM* x {t) - CiNxit). (376) 

Inserting (36) into (29), carrying out the integration and taking s = and s = — 1, the 
W-function and Q-function for the single-mode case can be obtained: 



WW(a u s,t) 



2(n\m\) exp 



e(t) ie(t) , 



e(t) ie(t) , 



(aie 2 -aje 2 ) 2 ( aie ~ 2~ +a*e — 2~ ) 2 



2(r(t)-s-2| 1 /.(t)|) 2(r(t)-s+2|V(t)|) 



TrVCrW-s) 2 -^*)! 1 



n h n I2 m+l2—n m+l\—n 

x E E E E E E (-2)—+-+-^ 



I 1= 0ni=01 2 =0mi=0 fei=0 fc 2 =0 

r!(m — n + mi + / 2 — A4 — r)![ni!mi!(ra — n + / 2 — fci)!(m — n + ^ — k 2 )\] 



h\k 2 \(h - m)!(/ 2 - mi)!(n - /i)!(n - fe)! 
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*[<(*) -<(*)]^^[«*) -«)]^nci(*) -ci(t)]*[C(t) -c(*)f 



xh-WCW + i7 + (*)C(*)] n " 1 [<(*)C-(*) + <(*)C + (*)1 



n-lo 



X 



ry;(t)c;(t) -<(t)C(t) r( y+(t)C + (t) -y-(t)Ut) 
{ z(t) ) l z(t) 



n\—r 



x[l-2(|C + (t)| 2 + |C-W| 2 )] 



2m- 2n+ij +^2 ~~ ^1 — ^2 



nj +m;L 



x[l-2(Mt)|^ + |»7_(t)| 2 )]^ 



x[i-^ic-(^-(*)-c + (^ + wi 2 r 



x Hi n 



xHi 



X(t) 



.y/ftit) -rji.it) 

Y(t) 



"l2—m\ 



X*(t) 



l WCl(t)-C(t) 



H, 



k-2 



Y*(t) 



y/cm - a 2 (t) 



. p(|mi— ni|,|ni+n+fci— m— Z2I) 



where 



^ 2 (t)-4|C-(t)r7-(^)-C + (% + (t)p 

^(t) + 4|C-(*)»/-(«)-C + (*)^(*)l s 



|r(t)| 2 = li^t)! 2 + iL^t)! 2 + iM^t)! 2 + lA^t)! 2 , 



i!>(t) = #i(t)Li(t) + N 1 (t)M 1 (t) = \^{t)\e 



ie(t) 



v/2(r(t)-s±2|^)|) 



u.(t) 



<(') 



C±(«) 



X(t) 



M*(t)e 2 ±N 1 (t)e 

v/2(r(t)-s±2g(t)|) ' 
r7 + (r)(aie 2 + a x e 2 j ?7_(rj(a 1 e 2 — a x e 2 j 

V2(r(t)-« + 2g(t)i) + ^/2 = (r(t)- S -2|^(t)r' 



(38) 

(39a) 
(396) 

(39c) 
(39rf) 

(39e) 



r 



20 

y / £ x = C+(t)(«ie 2 + a*e 2 ) (_(t)( ai e 2 - g*e 2 ) 
y 2 (r(t) - s + gggg y/2(r(t) - 3 - 2\j,{t)\) ' 

*(*) = >/[l " 2 (M*)I 2 + to-(*)l 2 )][l " 2(|C + (*)I 2 + IC-WI 2 )], (39^) 

= \[n\ + mi — 1 7&i — m-il], i? m is the Hermite polynomial of order m and P r (#) is the 



Jacobi polynomial which is defined as 



^w=D-i) ( -"(;!3( 

j._n \ / \ 



r + k + c + d\ ,x + 1,, 

k . (40) 



k r 2 ' 

Equation (38) is real in spite of its complex form, which can be seen explicitly in the 
summations where we can find each term with its complex conjugate. 

We can easily check the limits of equations (36) and (38) as t — ¥ 0, which give the 
corresponding well-known quantities for the Fock state \n) appropriate for the description 
before the interaction starts. In fact, this is clear also from the solutions of the Heisenberg 
equations of motion, where at t — the all factors reduce to zero except K\ (0) which equals 
1. So we get 

C( 1 )(C, S ) = exp[i( S -|Cr)]L n (|Cr), (41) 

W(a) = -(-l) n exp(-2H 2 )L„(4H 2 ), (42) 

71 

1 \a\ 2n 
Q(«) = -^exp(-|a| 2 ), (43) 

which are the s-parametrized characteristic function, W-function and Q-function for the 

Fock state \n). 

As we have mentioned before, the nonlinear directional coupler can be used as a source 



of quantum states 48j. This may be illustrated by displaying one of the quasiprobability 



functions J49| . The best quasiprobability functions for this task are W- and Q-functions since 
they are not singular and may contain oscillatory fringes (particularly VF-function) that are 
indicative of nonclassical behaviours. So we have plotted W-function and Q-function using 
(38) in Figs. 7 and 8, respectively, against x = Recti and y = Im«i, when the first state 
is the Fock state |1) and the second one is the vacuum state |0), i.e. n — l,ra = 0; A3 = 
1, Ai = A 2 = A 4 = 0.25 and for shown values of time. We have considered quasiprobability 
functions at t 6 [0,7r]. In Fig. 7a we have the W-function for t = -^, i.e. after short time 
interaction between the two modes we observe similar behaviour as for the W-function of the 
state |1) (see Fig. 2 of [50j]), which means that pronounced negative values are exhibited. 
This behaviour of the W-function is completely different by increasing the time (t = |); 
we see disappearance of negative values of the quasidistribution and a stretched positive 
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( 7a) 




( 7b) 




(a) 



(b) 



: 7c) 




(c) 

FIG. 7: W- function for the single mode (mode 1) for different values of time t when both the 
modes are initially in the number states; the first mode is in the state |1) and the second mode is 
in the state |0) and A3 = 1, Ai = A2 = A4 = 0.25: a) for t = 7r/100; b) for t = tt/2; c) for t = ir. 



iieak occurs (Fig. 7b). This form of W- function is close to that of squeezed vacuum states 



171 ]. i.e. squeezed vacuum states can be generated, in principle, in our model. It should 



be borne in mind that the specific direction of stretching for the quasiprobability function 
of squeezed states may be achieved by choosing a suitable value for the phase of squeezing 
parameter. Of course, in Fig. 7b, there is a negligible spike at the top of the peak which can 
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be smoothed out by governing the coupler parameters. After larger time interaction t = n, 
the negative values are reached again but they are less pronounced and asymmetry can be 
observed due to stretching (Fig. 7c). So we meet a time development of the W- function as 
a result of the power transfer between the two modes inside the coupler. The Q-function is 



( 8a) 



( 8b) 




(a) 



(b) 



( 8c) 




(c) 

FIG. 8: Q-function for the single mode (mode 1) for different values of time t when both the 
modes are initially in the number states; the values of parameters are as in Fig. 7. 



the quasiprobability function which is always positive definite, however, it can be used as 
an indicator for the squeezing in the model by including stretching in the phase space. In 
Fig. 8 we can see a kind of relation of the behaviour of W-function and Q-function and we 
can observe the top hole peak for short and long time interaction, which does not appear 
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for intermediate interaction times. For all cases the stretching is remarkable. 
(ii) Input coherent light 

In a similar way as we followed in the case (i) we can study the same quantities when both 
the modes are initially in coherent states. In this case the density operator is given by 
Pcoh(O) = \u)i\a) 2 2{a\i{aY (44) 

Then the two-mode s-parametrized characteristic function is derived in the form 



cS(Ci,C2,M)=exp{X; 



j = l L 



-5 



+ 2|L,(t)| 2 + 2|iV,(t)| 2 ) 



I0I 2 } 



x 



ex p{E 2 K www + mwM + c - c -] 



x exp CiC2[K{{t)N*{t) + M*(t)L* 2 (t)} 



c.c. 



x exp{-CrC2[^(*)Wi(*) + J\£(t)Li(*)] - c.c.} 



X 



exp{^ [&«*(*) -(*«(*)]}, (45) 

where ctj(t) are the mean values of the operators a>j(t) with respect to the coherent states. 
Therefore the two-mode s-parametrized quasiprobability function is 
W^( ai ,a 2 ,s,t) = 1 [(llaWf + lNtf)] 2 - |i?i(t)| 2 ) (5 + (t)5_(*) - ?*(«))]"* 



x exp 



x exp 



ff-^X 2 (t) + g+(t)X^(t) + 2X-(t)X+(t)T(ty 
S + (t)S-(t)-T*(t) 

^[El(t) + £* 2 (t)] - [|Li(*)| 2 + I^WI 2 ]!^^)! 2 



l^xC*)! 2 + |^i(*)| 2 - l^i(«)l 2 

where we have used the following abbreviations 
Mt) = ±(l-s + 2\L J (t)\ 2 + 2\N J (t)\ 2 ), 
Bj(t) = N*(t)M*(t) + L*(t)K*(t) = \Bj{t)\e 2i5 ^\ 
D(t) = K*(t)N*(t) + M*(t)L*(t) = |Z>(t)|e*<*>, 
C(t) = L 2 (t)iVi(i) + NZWlnd) = \C(t)\e*M, 
E 3 (t) = {^(t) - a 3 )e^\ 

F ± (t) = D(t) shift (t) + S 2 (t) - X (t)] ± C(t) sinft(t) - S 2 (t) + 7 (t)], 
i? ± (t) = £>(*) cosft(*) + <J 2 (t) - x(*)] ± C-(t) cosft^) - <J 2 (t) + 7 (t)], 



(46) 
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s4t) = Mt)-\B 2 (t)\ R - {t) F2{t) 



A 1 (t)-\B 1 (t)\ A^ + lSiWI' 
R-(t)F + (t) R + (t)F-(t) 

U Ax^-IB^)! A 1 (t) + \B 1 (t)\' 

x m - i\F m + P*mi + ^(flffifl-flift)] - MMgl±M)J 
x + (t) - »[^,(t) + E 2 (t)} + Ai(t) _| Bi(t) | - *— a^tJb^jT' 

X_(t)-[E 2 (t)-E 2 (f)} + Mt) _ lBM +* Al{t) + \ Bl (t)\ ' (47) 

with the following condition |A,(£)| > |Pj(t)| for the Glauber P-function and no additional 
constrains. 

From equation (46) we can see that W^(ai, a 2l t, s) includes the nonclassical correlation 
nature due to the presence of the terms aia 2 ,a{a 2 , etc. These mode correlations have 
been used in a number of studies on nonclassical aspects of light including questions like 



violations of Bell inequalities 5l|. The amount of correlation between the waveguides inside 
the coupler is governed by the coupler parameters, i.e. ctj,\j,t. Further, the P-function 
does not exist for |A,(£)| < \Bj(t)\, and this should be reflected as a nonclassical effect in 
the behaviour of the compound modes inside the coupler. The physical reason for this is 
that the modes may no longer fluctuate independently in even small amount allowed in a 
pure state. 



parametrized quasiprobability function are given, respectively, as 

-Ul- 8 + 2|L 1 (t)| 2 + 2|iV 1 (t)| 2 ) Id 2 + CaJ(t) - C«i(t) 



For the single-mode case the s-parametrized characteristic function and the s- 
•ametrized c 
C«(C,M) = exp 

x exp S^ewmMm + lk^w] 

x exp !.^C 2 [Ni(t)M 1 (t) + IatyKxit)] } , (48) 



^[V + l^iWI' + lMWPP-lBxCt)!' 
[1^ + \ Ll (t)\ 2 + iN^mia, - ai \ 2 
X<X1>1 r ^ + |L 1 (t)|2 + |AT 1 (t)|2]2_ | Bl (t)|2 

v f WB^WElJt) + Ef{t)] \ 

eXP l [^+\L 1 (t)\ 2 + Wi(tW-\B 1 (t)\^' l4yj 

and |Li(t)| 2 + |iVi(£)| 2 > |Pi(£)| must hold for the Glauber P-function. It is known that 

the correspondence between quantum and classical theories can be established via the 

Glauber-Sudarshan P-representation. But the P-representation does not possess all the 
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properties of a classical distribution function for quantum fields. More precisely, light 
fields for which the P-representation is not a well-behaved distribution (in most processes 
in interaction at least for some values of interaction time, including the process under 
consideration) exhibit nonclassical features. From (49), the P-function is not well defined as 
an ordinary function for |Li(t)| 2 + |iVi(£)| 2 < |-Bi(£)| and hence the nonclassical effects, e.g. 
squeezing of vacuum fluctuations and sub-Poissonian statistics can occur, as we have seen 
before. Furthermore, the nonclassical effect, especially squeezing of vacuum fluctuations 



( 9a) 







3 » 



( 9b) 




(a) (b) 

FIG. 9: W-function for the single mode (mode 1) for different values of time t when both the 



modes are initially in the coherent states; |oi| 2 
for t = 7r; b) for t = 2ir. 



02 



2 and Afc are the same as in Fig. 7: a) 



in the case of our system, can be recognized in the behaviour of VF-function (and/or 
Q-function) in phase space as shown in Fig. 9 for shown values of parameters. For t — 0, 
i.e. when there is no interaction between the two modes, the W-function is identical with 
that shown for a single mode representing a symmetric Gaussian bell in phase space. As 
soon as the interaction switches on (t > 0), we observe that the Gaussian centre is shifted 
and the rotationally symmetric function of the initial state at t — gets to be squeezed 
in various phase space directions in dependence on time, as demonstrated in Figs. 9a,b. 
In other words, the initial symmetric contour of the VF-function has been stretched as 
the interaction switches on, i.e. noise ellipse characterizing squeezed light appears, which 
rotates in the phase space as the interaction time progresses. 
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(iii) Input thermal light 

Signal beams are usually accompanied by thermal noise, so that examination of quantum 
fields with thermal noise is an important problem from both theoretical and practical points 
of view. Such thermal field can be generated by a thermal source composed of many in- 
dependent atomic radiators and consists of the superposition of waves of many different 
frequencies within some continuous range. These waves can be regarded as independent 
waves with random phases [16J. This field possesses uniform phase distribution (it is de- 
scribed by normal distribution), exhibits thermal statistics, i.e. g^(0) = 2, and its photon 
distribution is the Bose-Einstein distribution. Here we study the quasiprobability functions 
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FIG. 10: The cut of the VF-function for the single mode (mode 1) when both the modes are 
initially in the thermal states with n\ = fi2 = \/2, Xi = A2 = A4 = 0.25 , A3 = 1 and £ = §• 



for two modes as well as for a single mode as before, when both the modes are thermal. In 
this case the density operator takes the form 

fr(0) = 7Z -TTTZ — T J2 MXT MXT WlHwHlR (50) 

(ni + l)(n 2 + l) n ^ \ni + lj \n 2 + lj 
where fi\ {n-z) is the average thermal photon number for mode 1 (2). It is clear that the 

thermal distribution has a diagonal expansion in terms of the Fock states. This diagonality 

causes the electric field expectation value inside the coupler to vanish in thermal equilibrium 

at all times. This, of course, is related with the linearity of the relations (5) in terms of 

creation and annihilation operators. 
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The two-mode s-parametrized characteristic function is given as 
C t ( ?(Ci,C2,s,t) = exp 



"(ni + |)l^i(^)| 2 " {fh + ±)|Z 2 (t)| 2 + |(|Ci| 2 + |C 2 | 2 ) 



where 

Z 1 (t) = CiKt(t) - CiLi(t) + ( 2 M*(t) - C 2 N 2 (t), 

Z 2 (t) = ( 2 K* 2 (t) - QL 2 (t) + CiM*(t) - CrJVi(t). 
Therefore the two-mode s-parametrized quasiprobability function equals 



(51) 

(52a) 
(526) 



7rV[^(*)-|Ci(*)l 2 ][^(*)-|C 2 (*)| : 



x cxp 



x exp 



d(t)al + C{(t)al 



\D 1 {t)-a 2 \ 



2[A\{t)-\C 1 {t)\ 2 \ [^(*)-|C 2 (*)| 2 ] 

Q(t)[Di(<) - a 2 ] 2 + C 2 (t)[Dl(t) - a* 2 ] 2 



where we have defined 



2[%(t) - \C 2 (t)\*] 



(53) 



Mt) = (fk + ^)pi(t)| 2 + |^(*)| 2 ] + (na + i)[|Mx(t)| 2 + liV^t)! 2 ] - I, 

Mt) = (ni + ^)[|M 2 (t)| 2 + |iV 2 (t)| 2 ] + (n 2 + i)[|K 2 (t)| 2 + |L 2 (t)| 2 ] - i, 
(fk + l -)L\(t)Km + (n 2 + ^)M;(t)N*(t) 



Ci(t) 

c 2 (t) 



1 



A?(t)-|Ci(*)l s 



[C^lfit) + C*(t)l 2 2 (t) - 2A 1 (t)l* 1 (t)l 2 (t)] 



+2 



(n 2 + l -)L 2 {t)K 2 {t) + (n x + ^)M 2 (t)N 2 (t) 



D 1 (t) 



1 



{.ii(t)[ai/i(t) + all* 2 (t)\ ~ ail* 2 (t)Ci(t) - alk(t)Cl(t)} , (54) 



A?(t)-|Ci(t)l 2 
such that |Aj(t)| > |C,-(t)|. 

In equation (54) we have defined l\(t), and / 2 (t) as follows: 

h(t) = (ft! + n 2 + l)[Ll(t)M*(t) + KZ(t)NZ(t)], 

l 2 (t) = (m + l)[Ll(t)N;(t) + K*(t)M;(t)} + (n 2 + \)[Ml(t)K 2 {t) + Nl(t)L 2 (t)]. (55) 
We can see from (53) that the thermal light (classical light) propagating through the system 
under consideration can exhibit nonclassical effects, since the P-function can be singular 
under some constrains. Further we can see also that the nonclassical correlation between 
modes is available. 

For the single-mode case the s-parametrized characteristic and quasiprobability functions 
are 

.(1), 



C^(C,s,i) = exp 



-\C\V(t) 



2 ) ^ 2 



c 



u*(t) 



(56) 



then the s-parametrized distribution function can be written in the form 
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-(I) 



Kh(a,s,t) 



rvW)-§J 2 -|£W 

x exp 1 [J(t)-l?-\u(tw i ' (57) 

where we have denoted 

U(t) = L* 1 (t)K* 1 (t)(2n 1 + 1) + M*(t)N*(t)(2n 2 + 1), 

J(t) = [|Li(i)| 2 + l^iWI 2 ]^ + [iM^t)! 2 + 1^(01% + \ + |Li(£)| 2 + liV^t)! 2 , (58) 

with[J(t)-f] 2 >|f/(t)| 2 . 

It is well-known for the thermal optical cavity that photons have tendency to bunch 
each other, when photon distribution is described by the Bose-Einstein distribution (super- 
Poissonian statistics). However, as we have shown in section 3 the single mode thermal light 
can display squeezing of thermal fluctuations under this interaction, e.g. one can derive 
that the coupler is the source for squeezed thermal light. This can also be recognized in the 
behaviour of W^-function (see Fig. 10 where the cut through the jy-function is displayed). 
In this figure one can see the noise ellipse for squeezed thermal light with the center at the 
origin. 
6. Conclusions 

In this paper we have examined the quantum statistical properties of radiation gener- 
ated and propagated in the nonlinear optical coupler composed of two nonlinear waveguides 
operating by the second subharmonic processes, coupled linearly by evanescent waves and 
nonlinearly by nondegenerate optical parametric process. We have demonstrated regimes 
for generation and propagation of nonclassical light exhibited by squeeezing of vacuum fluc- 
tuations and/or antibunching of photons (sub-Poissonian photon statistics). We have also 
obtained quasidistribution functions for the initial light beams which are in coherent states, 
Fock states and thermal states. Compared to earlier results for nonlinear optical couplers 
we have shown that the nonlinear coupling increases in general quantum noise in the device 
even if in some cases it can support generation of nonclassical light. 

The motivation for examination of the system under consideration arises from the pre- 
vious investigations of the nonlinear couplers as promising devices to produce nonclassical 
light. When coherent light is injected initially in the system, squeezed as well as sub- 
Poissonian light can be generated. For injected number states, squeezed vacuum states are 
produced. When thermal light initially enters the coupler, the coupler can operate as a 
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microwave Josephson-junction parametric amplifier [37|. These effects have been recognized 
to result from the competition between linear and nonlinear properties of the system and 
are dependent on the initial amplitudes of the input fields. The crucial role plays here the 
mechanism of the energy exchange between waveguids. 
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